arXiv:1509.07886vl [hep-th] 21 Sep 2015 


Systematic Renormalization of the 
Effective Theory of Large Scale Structure 


Ali Akbar Abolhasani,“’^ Mehrdad Mirbabayi,'^ and Enrico Pajer'^ 

“ Physics Department, Sharif University of Technology, Tehran, Iran 

^School of Physics, Institute for Research in Fundamental Sciences (IPM), Tehran, Iran 

^ School of Natural Sciences, Institute for Advanced Study, Princeton, NJ 08540, USA 

Institute for Theoretical Physics and Center for Extreme Matter and Emergent Phenomena, 
Utrecht University, Leuvenlaan 4, 3584 CE Utrecht, The Netherlands 


Abstract 

A perturbative description of Large Scale Structure is a cornerstone of our understanding of the 
observed distribution of matter in the universe. Renormalization is an essential and defining step 
to make this description physical and predictive. Here we introduce a systematic renormalization 
procedure, which neatly associates counterterms to the UV-sensitive diagrams order by order, as 
it is commonly done in quantum field theory. As a concrete example, we renormalize the one-loop 
power spectrum and bispectrum of both density and velocity. In addition, we present a series of 
results that are valid to all orders in perturbation theory. First, we show that while systematic 
renormalization requires temporally non-local counterterms, in practice one can use an equivalent 
basis made of local operators. We give an explicit prescription to generate all counterterms allowed 
by the symmetries. Second, we present a formal proof of the well-known general argument that 
the contribution of short distance perturbations to large scale density contrast 6 and momentum 
density 7r(fc) scale as k'^ and k, respectively. Third, we demonstrate that the common practice of 
introducing counterterms only in the Euler equation when one is interested in correlators of 5 is 
indeed valid to all orders. 
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1 Introduction 


A robust and accurate understanding of the gravitational clustering of Dark Matter is one of the 
main goals of cosmology. The smallness of the smoothed density and velocity on large scales sug¬ 
gests that some appropriate perturbation theory should converge to the right answer in this regime. 
It has recently become clear [1, 2, 3] that a consistent treatment of this problem requires the addi¬ 
tion of an infinite hierarchy of effective interactions to the fluid equations of Standard Perturbation 
Theory (SPT) [4].^ Although the lowest order effective interactions have been discussed in details 
(in the density [2, 3, 17, 18, 19] and velocity [20] power spectrum, and in the density bispectrum 
[21, 22], including non-Gaussian initial conditions [23]), a complete prescription to generate all the 
relevant operators to all orders in perturbation theory is still missing. In addition, the current 
renormalization procedure - a procedure that makes the perturbative expansion physically mean¬ 
ingful - has been carried out only at the lowest order for a limited set of interesting observables, 
and not in a fully systematic way. 

In this paper, we introduce a renormalization procedure that works systematically to all orders 
in perturbation theory and give an explicit prescription to construct all effective counterterms 
allowed by the symmetries. Here we summarize our main findings: 

• Systematic renormalization In SPT one solves the continuity and Euler equations order 
by order: 

5 + e = - di{5v^), 

0 + ne + = - diiv^djv^). 

The nonlinear terms relate short and long wavelength modes to each other and lead to “loop” 
diagrams in the perturbative calculation of correlation functions. To be meaningful, loops 
need to be regularized and new interactions or “counterterms” must be included. With the 
exception of [20], all treatments so far have considered counterterms only in the Euler equa¬ 
tion. Moreover, the SPT kernels (and the associated diagrams) combine several qualitatively 
different contributions together, making it practically impossible to separate embeddings of 
lower order loop diagrams into a higher order diagram. Because of this, the connection among 
various terms in the loop expansion, e.g. 1-loop 2-point correlator and 1-loop 3-point correla¬ 
tor, is less clear than in the common treatment of for example quantum field theories. Here, 
we introduce a systematic renormalization that avoid these drawbacks and is analogous to 
the one commonly adopted in quantum filed theories (including the Peynman diagramatics). 
For additional clarity, the procedure is illustrated using the explicit examples of 1-loop power 
spectrum and 1-loop bispectrum. 

Our systematic renormalization makes explicit the implications of the long memory of the 
short-scale modes, of the equivalence principle and neatly organizes the generation of new 
counterterms. Nevertheless, it has its own drawbacks. First, it requires the introduction 
of counterterms in the continuity equation. These are indeed necessary to discuss velocity 
correlators, but they are degenerate with those in the Euler equation if one is only interested 
in correlators of <5. Second, systematic renormalization requires that all counterterms are 
non-local in time. A practical remedy to these shortcomings is presented in section 6. 

• Mass-weighted velocity/momentum description To renormalize the one and two-loop 
power spectrum and one-loop bispectrum, it is sufficient to include counterterms to the Euler 

^The discussion has by now been extended beyond dark matter [5, 6, 7, 8, 9]. For a recent discussion of the range 
of validity of the effective theory see [10, 11, 12). See also [13, 14] and [15, 16] for a different but related approach. 
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equation [2, 20, 18, 21, 22]. It is well-known that this approach works to all orders and for all 
statistics of 5, and is equivalent to a particular definition of velocity in terms of momentum 
current.^ In section 6.1 we summarize and sharpen those arguments. 


• All-order counterterms We provide an explicit prescription to generate the connterterms 
allowed by symmetries to all orders. The discussion separates into the momentum description 
and the systematic renormalization approach. The latter case is more cumbersome, so we 
refer the interested reader directly to section 4.4. The result in the former case can instead 
be stated concisely. The Euler equation must be supplemented by 


1 

TT6 




where the two-index symmetric tensor t ^ is constructed order by order by traces and products 

~ (n) 

of locally observable n-th order, two-index tensors and their spatial derivatives. These 
tensors (also called the “Eulerian” basis) are defined recursively by [8] 


= didj4>{x,T), and = 


(n) _ 




(n — 1)! 


(nj; '))'-(n-i)n 


(n-l) 


where prime denotes the convective derivative (4.10) with respect to logH {D is the linear 
growth factor). For example, at lowest order in derivatives and using matrix notation, the 
first few terms are 


E* lTr[n(^)],fl(^) 

2"'^ lTr[(n(^))2], l(Tr[n(^)])2, n^^^Tr^^^)], (n^^^)^ 

lTr[(n(^))3], lTr[(n(i))2]Tr[n(^)], l(Tr[n(i)])3, ITr^^^^fl^^)],... 


Two comments are in order. First, starting from quartic order, these operators can not be 
written as local functions of (second derivatives of) the Newtonian and velocity potentials 
(j) and 4>y [8]. Second, after taking the two derivatives of appearing in (1.1), some terms 
become degenerate. For example, there are five second order terms above but only 3 inde¬ 
pendent counterterms at this order [21, 22]. 

• Double softness It is well-known (see e.g. [24]) that the conservation of momentum and 
locality implies that interactions among modes of characteristic wave-number q can generate 
a longer wavelength mode k q suppressed at least as 6{k) oc k'^, namely the density is 
“double soft”. Analogously, the momentum is single soft, [tt] oc k. After reviewing a general 
argument for this property due to Peebles, we prove it directly from the equations of motion 
in section 5 (details are in appendix B). Double softness can be seen and implemented most 
straghtforwardly using mass-weighted velocity or momentum. On the other hand, in terms of 
the more general velocity used in [20] and in the systematic renormalization, double softness 
implies a highly non-trivial relation among counterterms to all orders. 

In this paper, we have extensively used the same nomenclature as in quantum and statistical 
field theory to stress the similarity with the problem at hand and to make the discussion more 
intuitive. When needed, we have adapted and refined some of the dehnitions. For the convenience 
of the reader, a glossary of the technical terms is provided in appendix E. 

^It should be stressed that care is needed when comparing the prediction for velocity correlators with simulation 
or observations. We refer the reader to [20] for a discussion of this issue. 
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The rest of the paper is organized as follows. In section 2, we review standard perturbation 
theory in a compact “vector” form and introduce diagrammatic rules, which we elucidate for the 
case of the tree-level bispectrum. In section 3, we discuss how interactions give rise to loops and 
the separation between 1-particle irreducible (IPI) and 1-particle reducible (IPR) diagrams and 
stochastic and non-stochastic counterterms. The one-loop power spectrum is used as an example. 
In section 4, we present our systematic renormalization procedure and discuss the implication of the 
long memory of the short-scale modes and of the equivalence principle. The one-loop bispectrum is 
used as a concrete example of our general remarks. Section 5 is devoted to the discussion of double 
softness. Finally, in section 6, we give the explicit form of all allowed counterterms and put the 
density-only description on a firmer footing. Some technical details are left to the appendices. 

2 Preliminaries 

At very large scales dark matter can be treated effectively as a pressure-less fluid of density contrast 
6 = {p — p)/p, and velocity u*, whose dynamics (in a matter dominated model of expansion rate 
%) are governed by the continuity and Euler equations: 

Ue = -di{5v^), 

V + T-iv + V<j) = — V ■ Vv, 

where the gravitational potential satisfies V'^cj) = Starting from zero vorticity, i.e. V x u = 0, 

it remains zero under this evolution. Hence, we can take the divergence of the latter equation and 
use the velocity divergence as the independent variable 0 = SjU*: 

e + W+ ^71^5 = -diiv^djv^). ( 2 . 2 ) 

(We will comment on the generation of vorticity in section 4.3.) At shorter scales the above 
description is incomplete and higher derivative corrections must be included. In practice, the need 
for such corrections (counterterms) can be seen from the dependence of the nonlinear solutions on 
the choice of the ultraviolet cutoff. The coefficients of the counterterms depend on the cutoff in such 
a way that the final result is cutoff-independent, though it does depend on the unknown details 
of the ultraviolet physics. That is, the counterterms capture the contribution of short distance 
physics by introducing new effective interactions among long-wavelength modes. Their coefficients 
are to be fixed either by a first-principle calculation or, more realistically, using simulations and 
observations. In this paper, we will derive the general form of these counterterms. 

For this purpose, it is useful to develop a systematic approach to renormalization—the procedure 
of removing cutoff dependence—as in conventional relativistic field theories. The main difference 
is that here we are dealing with an initial value problem, whose implications becomes clear in the 
sequel. In particular, it helps to have a more symmetric treatment of the primary variables <5 and 6, 
and a more explicit diagrammatic representation of perturbation theory. Throughout the paper we 
work in an Einstein-de Sitter (EdS) universe but to a good approximation the conclusions remain 
unchagned in ACDM cosmology. Following [25, 26] we introduce a doublet field 

The equations of motion can then be written (in momentum space) as 
^ /* 

—'^a{k,r]) + il.ab{v)^b{k,rj)= 7^^),(/c, fci, ^ 2 ) ^^(fci, ?/)^c(fc2, ??), (2.4) 

JkiM 
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where 


r 

r d?k 


0 -1 

Ik 

J (2vr)3 ’ 

77 = log a, Qab = 

-3/2 1/2 


and non-vanishing components of the symmetrized vertex functions 7 ^®^ are given by 

7i2i(fc, fci, A:2) = - fci - , 

7ii2(fc, ^ 1 ,^ 2 ) = 5^{k-ki- , 

7222(fc, ki,k2) = 5^{k - ki - fc2)/3(fci, ^ 2 ) • 

The vertices can be read from (2.1): 


a{ki,k2) 


ki ■ {ki + k2 


I3{ki,k2) 


ki ■ k 2 (fci -I- ^2)^ 

2 kfk2 


(2.5) 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 


2.1 Perturbation theory and diagrammatic rnles 

Given some initial condition ^a{k,'r]o) = {5i{k),6i{k)), the above system can be solved perturba- 
tively to find 'I'a(/c,T). At first order, we have 

^ a \ k , r ]) = T)+(77,770)^6(^,770), (2.10) 

where repeated indices are summed over and is the growing solution of the linearized system. 
Note that it is fc-independent. D^^{ri) is related to the conventional growth factor through 


D(v) = Vo) = Y ^ 2 aiv, Vo)- (2-11) 

a a 

The higher order solutions can be expressed, recursively, by combining two lower order ones via 
the interactions vertices 'jabc- 


= Y [ [ dv'vabiv-v')'ybcd{k,ki,k2,r]')^lrHf^i^v')'^f ""^^2,^'), ( 2 - 12 ) 

m=l -Im 


where the retarded Green’s function, or propagator, is given by 


9ab{v,v') 


5 


- 1 

CO 

to 

g-fG-V) 

1 

1 _ 

2 

3 2 

' 5 

1 - 

CO 

- 1 

CO 

1 


(2.13) 


for 77 > 77 ', and 0 otherwise. The growth function Zl)))j( 77 , 770 ) coincides with the growing part of 
dabiViVo) and, since 770 is practically — 00 , it can be approximated by gabiv^Vo)- This perturbative 
expansion is diagrammatically illustrated by Feynman diagrams as 


gab{v,v')- — 

a, 77 


k 

- 

b, 77 ' 



(2.14) 
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Every propagator is represented by a line, and the cubic interactions by vertices. The arrows 
indicate the flow of time, and clearly the first index of 'jahc^ which corresponds to the out-going 
line, has no symmetry property with the last two which refer to the in-going lines. A typical 
diagram contributing to is 





(2.15) 


where here and in what follows the indices are often omitted to avoid clutter. Note also that 
the momentum integrals in ( 2 . 12 ) must be understood as running up to some ultraviolet cutoff A 
beyond which the effective theory is not applicable anymore. Renormalization ensures that the 
final results do not depend on the choice of A. 

Using the above Feynman rules, the following recipe can be used to calculate correlation func¬ 
tions perturbatively. An external leg is given by the sum over all tree diagrams constructed 
by successive bifurcations through trilinear interaction vertices to get to n initial fields 6i{k), as in 
(2.15). Any trilinear vertex merges a and a to result in an (n-l-m)-th order field. 

These perturbative expansions for the external fields are then correlated using the initial statistics 
of 6i{k). We assume Gaussian initial condition 

{6i{k)6i{k')) = PUk){27rf6^{k + k'), (2.16) 

which reduces the correlation functions of several 'h fields to a sum over all possible pairings of the 
initial fields that appear in their perturbative expansion. The correlation of initial fields is called 
contraction and is reperesented by a dot. If all initial momenta in a given diagram are fixed in 
terms of the external ones, the diagram is called tree-level. One example is 


(2.17) 


Otherwise, if there are undetermined initial momenta it is a loop diagram. There is one unde¬ 
termined momentum for each loop and it must be integrated over. The tree-level contribution 
to an N-point correlation function is of order 2{N — 1) in di, and an n-loop contribution is of 
order 2{N -|- n — 1). It is common to label various diagrams by the orders of their external legs 
in perturbation theory. For instance, the 1-loop contribution to the power spectrum where is 
correlated with is often called P 13 . 
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2.2 Example: tree-level bispectrum 


Let us conclude by illustrating the above formalism in the example of tree-level 3-point correlation 
function of density contrast 


{5{ki,r])5{k2,r])5{k3,r])) = {27rf5^{ki + k2 + k^)B{^2{^i^^2, k3,ri). 
In the conventional SPT formalism this is given in terms of the kernels: 

Bli 2 {ki,k 2 ,k 3 ,ri) = 2 e^^F 2 (fci, A: 2 )Piin(fei)^’iin(fc 2 ) + 2cycl.perm., 


where 


17 ( 1 , r. ^ ^ 2 (fci • fc2 

F2{ki,k2) ^ + ^ ^2^2 




, / 1 1 

^ 2 ( 7:2 + 1:2 


'kl 


kl 


In the doublet formalism, the relevant diagram is 


1 




(2.18) 


(2.19) 

( 2 . 20 ) 


( 2 . 21 ) 


where we demonstrated the SPT diagram on the right for comparison. Following the Feynman 
rules, one obtains 

X, / dv'giaiv, g) gbdiv, go) gceW, go) gifig, go) gigig, go)iabc{ki, k 2 , fc3)^’iin(^i)Piin(A;2)- 

abode f g ^0 


Using the above expressions for gab and jabc (and, of course, including other permutations), it is 
easy to show that this agrees with (2.19). 


3 Loops 

Loops are encountered in the perturbative calculation of correlators when the initial momenta are 
not fixed in terms of external momenta, and therefore, are integrated over. The ultraviolet part of 
these integrals include short wavelength modes which are not expected to be well described by the 
effective field theory. The integrals are often divergent. One should therefore regulate the integrals 
(say by sharply cutting them off at momentum A), and use the cutoff dependence to introduce new 
effective interactions (counterterms) that capture the contribution of the short scale physics to the 
low momentum observables as an expansion in k/A. This procedure provides a systematic way 
of inferring the form of higher derivative corrections to be added to the long-wavelength effective 
theory, and in most cases, it covers all possibilities that are compatible with the symmetries of the 
system under consideration.^ 

Hence, understanding better the properties of loop diagrams is the key to identify the general 
structure of higher derivative corrections. This is the subject of this section. 

®There can be exceptions when some higher derivative operator is protected by a non-renormalization theorem. 
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3.1 Reducible and irreducible diagrams 

Consider the P 13 contribution to the 2-point correlation function 



Suppressing time-integrations and indices the diagram would have the following momentum struc¬ 
ture 

Plinik) j (fqPiin{qhik,q,k - q)j{k - q,k,-q). (3.2) 

To renormalize this loop integral a linear counterterm must be added to the equations of motion. 
Its leading contribution to the 2-point correlation function 

—<— ® 4—•—►— (3.3) 

combines with the original diagram to make the full result cutoff-independent. 

Now suppose the same loop diagram of (3.1) is inserted into a higher order graph, e.g. 



(3.4) 


One would expect that if the counterterm (3.3) is chosen appropriately, inserting the same quadratic 
vertex in place of the loop diagram, namely 



should renormalize (3.4). We will see that such a choice of the counterterm is indeed possible, and 
in particular, the diagram (3.5) does renormalize diagram (3.4). Hence, embedding lower order 
loop diagrams in higher order ones does not lead to genuinely new corrections (counterterms). We 
call these loop diagrams reducible, or in analogy with particle physics diagrammatics, 1-Particle 
Reducible (IPR). Their characteristic property is that there exist an internal line I that is not 
immediately connected to any external line through a contraction such that cutting I divides the 
diagram into two disconnected pieces. All other diagrams are called 1-Particle Irreducible (IPI) 
diagrams, (3.1) is an example. 

Since the counterterms are to be added to the equations of motion as new vertices, it makes 
sense to work not with the full diagrams representing the correlation functions but directly with 
the time-evolution diagrams. These show the evolution of a set of initial fields into one higher order 
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final field at a later time. The loops correspond to integrating over some of the initial momenta, 
and hence, to allowing them to become “hard”, i.e. of the order of the ultraviolet cutoff A. We 
represent these hard momenta by thick lines. Moreover, to exclude IPR diagrams one would simply 
amputate all other (“soft”) lines. That is, the relevant diagrams are made of (i) several hard lines 
extending all the way to the initial condition (5i, (ii) one soft out-going line, and (iii) any number 
of soft in-going lines which are immediately connected to the hard lines. An example is 



(3.6) 


The presence of high-momentum hard lines makes these diagrams cutoff dependent. They are 
renormalized by adding new vertices to the effective theory that have the same structure of the soft 
lines, and the same cutoff dependence in their coefficients: 



(3.7) 


Note that the contribution of the new vertex to the final field always comes with the propagator 
of the out-going soft field which is common between (3.6) and (3.7). Since we are interested in the 
structure of the corrections to the equations of motion and not the resulting final field, we will also 
amputate the final propagator. 


3.2 Stochastic vs. non-stochastic corrections 

There is another important classification of counterterms that is unique to the case of initial value 
problem. The high-momentum initial fields were seen to result from integration over unfixed initial 
momenta. If in a given time-evolution diagram all of these can be paired and contracted with each 
other, the associated counterterm is called a non-stochastic one, e.g. 



The counterterms of these diagram capture the fact that short scales are affected by large scales. 
Otherwise, if at least one pair of the initial hard fields do not have exactly opposite momenta, it 
corresponds to the probability that short scale fluctuations average into a long-wavelength mode. 
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e.g. 


k 


stochastic diagram: 



(3.9) 


The associated counterterms are called stochastic. These are by nature indeterministic and hence 
contain a stochastic field (usually denoted by J{x,ri)) about which we can only make statistical 
statements. 

One can contract and integrate over the momenta of all initial hard fields that are paired in a 
time-evolution diagram (which include all of the hard lines in a non-stochastic diagram). This is 
because the loop is independent of how the ingoing soft lines are embedded in a more complicated 
diagram. For instance, the hard modes 6i{q) and Si{—q) in fig. 3.8 can be contracted regardless of 
whether we consider diagram of (3.1) or (3.4). After doing so these diagrams (and the associated 
counter-terms) incorporate the average response of the short modes to the initial long modes. ^ On 
the other hand, the unpaired initial hard modes will be contracted by unpaired modes coming from 
other external lines as in 


o 


(3.10) 


Our primary focus will be on non-stochastic counterterms since in our universe the stochastic ones 
are expected to be sub-dominant at long wavelengths (see below). 

3.3 Example: One-loop power spectrum 

As an explicit example let us consider the 1-loop non-stochastic contribution to the power spectrum 
(3.1). This corresponds to 


*-0 


◄- 7]', b 

'> 12 , 0.2 

whose amplitude is given by (recall that amputation removes the final propagator) 

rm 

/ dr ]2 02,112) ga 2 b{V 2 ,v') ■ 

Jri' 

Here 

Cairil-,k,b,ri 2 ) = 4c/a2 l(m, %)5fe2 ife, %)5ai fei (m, f/ 2 ) 

X / -Plin(g) Taaiaa +q, k - q) 'Jbi b 2 b {k - q,-q, k) , 

Jq 


(3.11) 


(3.12) 


(3.13) 


^Notice that the vertices (2.6) vanish for zero outgoing momentum since a{k, —k) = P{k, —k) — 0 (see section 4.3 
for more details). Therefore, non-stochastic diagrams with no ingoing soft lines all vanish. 
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the time variable is defined as a = e^, ga 2 i{vi^Vo) = fl'a 22 (??ij%) is the same as the growth factor, 
and the summation of repeated indices is assumed. Using the explicit expressions for the propagator 
and vertices, and restricting to q ^ k, we obtain the following expression for the ultraviolet part 
of the loop diagram 


Ca^{Vi-,k,b,r] 2 ) = e 

The velocity dispersion, which is defined as 



3 

-^9 

-4 

Sr?! 

6 2 + 


1 

5r?2 

6 2 


5 

5 


5 

5 

J 




= j dgPiin(g), 


(3.14) 


(3.15) 


is the cutoff dependent coefficient that we aim to renormalize by the addition of a counterterm 
proportional to k"^. This will be done in the next section. 


4 Renormalization and counterterms 


In this section, we hrst explain how to read off the counterterm associated to each amputated loop 
graph. Then, we discuss various properties of these counterterms, give their general structure, and 
conclude with an explicit application of these ideas to the 1-loop bispectrum. 

The amplitude of a generic amputated diagram (such as (3.6)) with the final vertex inserted 
at time r], outgoing momentum and index (fc,a), and a set of soft ingoing lines characterized by 
{kn,bn,r]n}, a^d hard ingoing lines characterized by {qi} (recall that hard lines are all initial fields 
5i{qi), so bi and ry' need not be specified), is of the form 


{dr]n}A 




[n;{kn,r]n};{qi}] 


'Vo 


'[{9a„bAr]n,r]n)'[{dl{q^). 


(4.1) 


In the non-stochastic case the hard lines {qi} can be paired and contracted to give Pimiqi) and 
integration over their momenta results in 


rv 

/ {dTln}Caai---a„ [V't {kn^Vn}] 9anbn ijlni Vn)- 
■^VO n 


(4.2) 


The Cab of last section was an example of this with only one ingoing soft line. The C kernel is a 
cutoff dependent quantity. This diagram is renormalized by the addition of a counterterm of the 
following form to the equation of motion (c.f. (3.7)) 


/ ' -I—r 

/ {dqn} C-aai---an [q; [K, qn}] || {K, qn)- 

''lo n 


(4.3) 


To cancel the cutoff dependence of the loop diagrams, the dependence of Ca... on {kn} has to match 
that of Ca.Y, namely the cutoff dependent part of Ca- when it is expanded in powers of k/q. Let us 
analyze this more closely. In the ultraviolet regime the momentum dependence of Ca— can always 
be brought via Taylor expansion into a sum of terms, each of which factorizes into a product of a 
function of {qi} times a function of {kn}. We are interested in the functional dependence on {kn}. 
Consider a vertex at qn where the soft line is attached to the diagram. From the equations of 
motion (2.1) and (2.2), the leading order vertices of the theory are 

— v^di6, —96, —v^diO, — {piV^) djv''. (4.4) 
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Thus, the soft line can realize one of the following six possibilities: 


^^2{kn,Vn), 


di6 ^ 9 ^ '^2(kn,Tjn)i 


5 ^ ^l{kn,r}n), di9 kn^2ikn,'nn), 


diV^ —7> 





(4.5) 


These are the only ways in which Ca... may depend on the field, up to possible multiplication 
by positive factors of k\, arising from vertices appearing at later times in the diagram. The 
contribution from these later vertices is guaranteed to be analytic in k"^ because (i) by definition 
of an amputated graph, a soft line n attaches immediately to a hard line, say with momentum q„, 
thus (ii) the later vertices depend on k^ only through + k^, and (iii) since -C Qn one can 
Taylor expand in fcjj. 

Moreover, Ca--- has to have a cutoff dependent coefficient with exactly the same time-dependence 
and overall amplitude as the one arising from the integration over {qj} in C^.Y■ However, there can 
be a cutoff independent piece with different time-dependence. This counterterm by construction 
renormalizes all possible embeddings of the original amputated diagram into any higher order 
diagram. Next, we will discuss the manifestation of characteristic properties of the cosmic fluid in 
the form of these new, and yet rather abstract, corrections. 


4.1 Long memory effect and non-locality in time 

Counterterms capture the effect of the short wavelength modes on the long wavelength physics. 
They result from replacing hard loop diagrams with new effective interactions among soft lines 
(i.e. integrating out the short-distance physics). The result of this procedure is guaranteed to 
be spatially local since the short modes can probe the long modes only in a region of size 
over which the long modes can be accurately Taylor expanded in powers of k/A. However, before 
reaching the virialization scale, the short modes evolve with approximately the same time scale 
as the long modes. Hence, we should expect the counterterms to be generically non-local in time 
as a short mode can be influenced by long wavelength modes at some earlier time and respond at 
a much later time. The hard lines in a loop diagram can therefore be collapsed to a point spatially 
but not temporally [15, 18]. 

This expectation is manifest in our general formula (4.3) and in the explicit 1-loop example. In 
order to systematically renormalize the amputated diagram of (3.12) one needs to add a temporally 
nonlocal quadratic vertex, corresponding to a linear counterterm in the equations of motion: 

rv 

dr^^aik,r])+nab^b{k,r]) = k'^ dqi Cafe(q;qi) Tb(fc,qi)H-, (4.6) 

Jo 

where we omitted the standard cubic vertices. Here k'^Cab has a cutoff dependent piece that (up to a 
minus sign) exactly matches in (3.14) so that the sum of the two cancels, up to a regularization 
independent piece. As will be seen in section 4.5.1 this counter-term automatically renormalizes 
the IPR contribution (3.4) to the bispectrum. On the other hand, if we were only interested in 
renormalizing the Pis diagram (3.1), this vertex could be replaced by a local one (the conventional 
speed of sound and viscosity counterterms) 

Ca(?7)A:^Ti(q), (4.7) 

where 

Ca(r/) =dvi Caa,{v,vi) D^bivi)- (4.8) 
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As we will see in section 6 this reduction is always possible in perturbation theory. However, it is 
incompatible with systematic renormalization. For example, using the local-in-time counterterm 
in (4.7), on finds that the IPR diagram (3.4) requires new counterterms [21]. 


4.2 Large scale flow and the Equivalence Principle 

The origianl system (2.1) is invariant under the addition of a uniform bulk flow [27, 28] 

x^x + nij), v^v + h{T), (4-9) 

where n is an arbitrary time-dependent vector. This ensures that a very long-wavelength pertur¬ 
bation, whose main effect on much shorter scales is to produce a uniform acceleration, is locally 
unobservable, in accord with the Equivalence Principle.^ Local observables can depend on the tidal 
field didj(j), the shear diVj, their spatial derivatives, and their convective time-derivatives 

Dr^dr + v-V. (4.10) 


The counterterms which characterize the response of the short scale physics to long-wavelength 
modes should also follow the same rules. The short modes can only respond to the locally mea¬ 
surable quantities. However, when the counterterms are non-local in time as in (4.6), they depend 
on these local observables along the large scale flow. The same considerations imply that the time 
integral must be taken along the fluid trajectory, which is implicitly defined by 

X{i[T']X,T] = X — j dT''v{Xfi[T"]X,T],T"). (4-11) 


In the following, we often drop the last two arguments of Xfi when no ambiguity arises. For 
concreteness, let us focus on the counterterm in (4.6). To be invariant under (4.9), it must be 
dressed to read [18] 

/ dr]' Cabimv) '^b{x{i{r]'),r]'). (4.12) 

The significance of this replacement x ^ Xf[ soon becomes clear in the example of 1-loop bispectrum 
in section 4.5.2. Perturbatively expanding x^ around x, we get 




rn 

,(£Cfi[r?';®,r?],r/') =Ta(r/') - / dr]iv{r]i) ■ V^aiv') 

Jri' 


+ 


/ drii / dr]2 [vim) ■ Vu(r/i)] • VTa(r/') -k 

Jr/' J rji 


(4.13) 


where all omitted space arguments are now x. This replacement essentially connects an infinite set 
of counterterms (and hence an infinite set of amputated diagrams) to one another. These are all 
diagrams in which the original diagram (3.8) is being displaced by the bulk motion produced by 

®For a discussion of these symmetries from the perspective of Weinberg adiabatic modes see appendix A of [21]. 
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any number of ingoing soft lines attached at later times than the original attachment of ki: 



(4.14) 


More generally, whenever a soft line realizes the part of the vertices in (4.4), it can be at¬ 
tributed to the shift from x to of a lower order diagram with the soft line removed. These 
diagrams have the following characteristic property in perturbation theory: In the presence of a 
long-wavelength mode with constant amplitude 5i{ki) and vanishing they diverge as k^/kf^, 
namely, proportionally to the bulk velocity induced by the long mode. This will be referred to as 
infra-red (IR) singularity. The UV part of these diagrams are fully fixed in terms of that of the 
original (unshifted) diagram. The counter-terms associated to these diagrams too are fully fixed in 
terms of that of the original diagram. Therefore, they do not introduce any new parameters into 
the effective theory. Eliminating these diagrams by shifting a; —>■ ajfj ensures that the new higher 
order counterterms needed to renormalize the remaining higher order IPI diagrams depend only 
on locally observable quantities, again measured along the dark matter flow x^. In other words 
there is no IR singularity, except in the argument Xf[ of the fields.® 

4.3 Soft outgoing momentum and vorticity 

Another property of loop diagrams to be discussed is their softness in the outgoing momentum. In 
perturbation theory, it can be easily verified that all vertices are at least single soft: 

\mi ’^abc[k,q,k - q) = 0{k/q). (4.15) 

q^oo 

In any amputated loop diagram, the latest vertex is one in which two hard lines combine to give 
an outgoing soft line of momentum k (see e.g. (3.6)). Because of (4.15) the UV part of this 
diagram will contain at least one factor of k, and so do the corresponding counterterms. Rotational 
invariance implies that the full amputated diagram must be of second order in soft momenta since 
there is no singularity (i.e. negative power) in the ingoing soft momenta after taking care of the 
shift terms as explained in the previous section. Thus expressed in terms of T(j, the counterterms 
must start from second order in derivatives, one of them being an overall derivative because of the 
property (4.15) of the final vertex. An example is 

( 9 i[ 4 'ai (®fl(r?i), r?i)Ta2 (®fl(?/2), ?/2)diTa3(®fi(r/3), r/3)]. ( 4 . 16 ) 

The ultraviolet part of Cabiq', k,r]), calculated above agrees with this expectation. 

The above softness property of the connter-terms has a natnral explanation based on the strnc- 
tnre of fluid equations. The continnity equation relates the change in mass density to the divergence 

®This fixed IR singular part can be used to derive a universal and dominant piece of the squeezed correlation 
functions, when one of the external momenta is much smaller than the others [29]. This piece often vanishes in 
equal-time correlation functions, however there is an exception in the presence of BAO feature [30]. 


15 




of the mass current (momentum density)^ 

5 = -V-TT, 


where tt = (1 + 5)v. 


(4.17) 


However, this relation between tt and v depends on the dehnition of velocity. Fixing velocity at 
some scale A to be given by [vr]A/[l + 5 ]a and then coarse-graining to a smaller cutoff A' results 
in a velocity field that differs from [ttJa'/)! + Hence, when working with velocity field as the 
primary variable one must add a locally observable vector as a counterterm to the mass current. 
The absence of IR singularity ensures that a locally observable vector must be at least of first 
order in derivatives (in the same sense that (4.16) is of second order in derivatives). Upon taking 
the divergence, this leads to a counterterm V • jc in the continuity equation that is second order 
in derivative counting, and because of the overall derivative it contains at least one factor of the 
outgoing momentum. 

Similarly, integrating out the short wavelength modes can introduce an effective force fc in the 
Euler equation. This is again a locally observable vector that incorporates the counterterms. After 
taking the divergence to obtain an equation for 9 it leads to a counterterm V • (/c/(l + S)) with an 
overall derivative. 

It is worth noting that since the original system of equations (2.1) do not generate any vorticity, 
to cancel the cutoff dependence of just SPT loops, fc /(1 + 5) can be taken to be curl-free. In other 
words, the counterterms to the 9 equation Ca=2,— only a total divergence but they are the 

Laplacian of locally observable scalars. We will verify this expectation in the explicit calculation of 
quadratic counterterms (see appendix A.2). Of course, the finite pieces of the counterterms do not 
have to obey this rule, and in principle the short scale dynamics can generate vorticity (see section 
6.1 and appendix B.3 for comments on the non-uniqueness of the definition of velocity). 

The property (4.15) and the fact that all counterterms have one overall derivative ensures that 
if the outgoing momentum k is much smaller than the soft ingoing momenta {ki} the solution for 
6{k) and 9{k) is going to be single-soft in k. In fact, as we will discuss in section 5, momentum 
conservation implies that 5 must not only be single-soft but rather double-soft in k in the limit 
k ki- As a result, only a subset of locally observable forces - those that respect this property - 
are allowed to be added as counterterms to the system. 

4.4 The list of counterterms 

Using the above rules we can write down an over-complete list of counterterms which can be added 
to the equations of motion. The softness in the outgoing line implies that they must be of the 
form did. d encodes the memory of all incoming soft lines that are attached to the hard lines of 
the amputated diagrams (3.6). Therefore, each field in d{x,r]) is integrated from rjo to r] against a 
general time-dependent memory function. Any attached ingoing line is either a locally measurable 
quantity made of the soft mode such as didjcj) or diVj (and their convective time derivatives (4.10)), 
or it is a velocity, which is not a local observable. As argued above the latter is necessary to 
perturbatively dress the time integrals at constant x to those along the fluid flow. In the former 
case, since both didjcj) and diVj have two indices, rotational invariance implies that to construct a 
vector d one needs to have one additional derivative. 

We proceed as follows. For simplicity, suppose there is zero vorticity; this allows us to write 
V = S/cpy, and introduce 

c/a = = {c/, c/y). (4.18) 

^Notice that we have changed the traditional definition of tt by a factor of 1/p, in such a way that tt is dimensionless. 
This is convenient because p completely disappears from all the equations. 
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Next, as an intermediate step we construct order two-index tensors made of various con¬ 
tractions of the local observables Ha = didj^a- Using matrix notation, the first few orders are 

uU n^iTi-n, (4.19) 

: n,Tr(nb), lTr(n,)Tr(nfe), lTr(n,nfe) (4.20) 

yd: n,Tr(nb)Tr(nfe),n,Tr(n,n6),n,nfeTr(ne),n,nfen„--- ( 4 . 21 ) 

where = 5®-^ , and as an example, Tr(nanfe) = Each 14(2 evaluated at a different 

time which is integrated against a temporal kernel until the final time. An over-complete (see 
below) list of vectors c® consists of all possible ways of acting by dj on any of Ila appearing in . 
Finally, one takes the divergence did. One second-order example is 

di / dr/idry2Araaia2(??,??i,f?2)9jn)2^^(aifl(r/i),r/i)n^2(®fl(^2),??2), (4.22) 

Jm 

where the derivative in djlia^ should, in principle, be with respect to aifl. However, the difference 
between this and the derivative with respect to x can be expressed in terms of higher order terms 
that are included in our list - in other words 

-—4- = (5®-^ -|- local observable tensors. (4.23) 

dx^ ^ ' 

So we take it to be djdx^ in what follows. Note also that the arbitrary time-kernel makes the use 
of convective time-derivatives redundant. 

Unfortunately the above procedure generates an over-complete list of counter-terms, as it allows 
us to introduce arbitrary observable vectors fc and jc as counter-terms. On the other hand, as 
we will argue in section 5, only a subset of these vectors - those which ensure double softness 
of 6 when kout "C {A:j}in “ are permitted by momentum conservation. An allowed list can be 
obtained by restricting the kernels to those which arise from perturbation theory. This elects 
certain combination of these vectors that are necessary to cancel the cutoff dependence of the 
loops. These combinations inherit the double softness property of the leading order perturbation 
theory which is discussed in section 5.2. However, in this way all amputated diagrams with fixed 
number of ingoing soft and hard lines combine to give only one new counter-term whose coefficient 
can be allowed to have a cutoff independent piece. While this is sufficient to cancel the cutoff 
dependence of the loops, in order to produce a complete basis that can give an effective description 
for all possible microscopic scenarios one needs to proceed to diagrams with increasingly large 
number of hard modes. 

We will next verify the adequacy of this formalism in cancelling the cutoff dependence in the 
explicit example of the 1-loop bispectrum. Several properties of the counterterms become apparent 
in this example. In section 6, we show how this rather cumbersome and nonlocal formulation can be 
reduced to a much more tractable and local one in any practical calculation by giving up systematic 
renormalization. Moreover, using the results of section 5, we will provide an alternative list of local 
counterterms that is complete but not over-complete. 

4.5 Example: One-loop bispectrum 

In this section we discuss the renormalization of the non-stochastic part of the 1-loop bispectrum 
which is conventionally called H 114 . The IPI and IPR contributions are evaluated and discussed 
separately. 
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4.5.1 IPR diagrams 

As expected from the general argument at the beginning of this section, it is easily seen that the 
IPR diagrams are automatically renormalized by the insertion of the linear counterterm in (4.6) in 
a tree-level diagram. The diagram 


ki,b-,ri 



(4.24) 


corresponds to taking the expectation value 

pTJ pTf\ 

/ dr^igaaA'n^Vi) / dr]2 Ca\a2 {m,V2;k) /4'i2)(A:,r?2)4'[^^(fci,r/)TW(fc2,??)) + 2 perm. (4.25) 

Jr}o Jr}o ' ' 

We denoted the UV part of this loop diagram by and it will be renormalized by adding the 
contribution of the diagram 


ki,b;r] 



(4.26) 


where the cross represents the counterterm in (4.6). The expression for this diagram is given by 
an identical expression as (4.25) with C —>■ A:^c, which have equal and opposite cutoff dependences. 
The other IPR diagram 
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has an expression 


/ rv rV2 rm 

(fq I drji / dq 2 / 5aai (^, m) ^? 2 ) 7ai bi ci(fc, ^ 2 ) C'fc2b3(r?2, ' 

Jvo Jm Jm 


(4.28) 


(fci, (fci, r?3)4'W (/C2, (^2, r/i)) + 5 


perm. 


Similarly, the UV part of this diagram is renormalized by the diagram containing the same 
counterterm in place of the loop: 


ki,b]rj 



(4.29) 


which has the same expression as (4.28) with C —>■ k'^c. 

It is worth emphasizing here that the temporal non-locality of counterterm Caa'iv]'n'^k) is 
crucial for it to cancel the cutoff dependences of all IPR loop integrals. This accommodates the 
possibility that the response of short modes should strongly depend on the time-dependence of the 
long modes. If only local in time counterterm are used, the UV divergences in IPR diagrams are 
not generically canceled by the lower order counterterms, spoiling the systematic character of our 
renormalization procedure. This is expected and in agreement with the general argument presented 
in section 4.1 about the long memory effect. 


4.5.2 IPI diagrams: shift terms 


As discussed in section 4.2, some of the IPI diagrams contributing to B 114 , whose (A:) field is 
of the form 



(4.30) 
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and 



(4.31) 


play the role of shifting the temporally nonlocal lower order counterterms to the fluid position (here 
2 indicates the velocity field). It is implicit in the above diagrams that the lines marked by 2 realize 
V m. V ■ V'I'a vertices, and that r ]2 > r/i. These correspond to contributions that diverge as l/k 2 
when k 2 —)• 0. We expect them to be renormalized by including the first shift term in (4.13) in the 
counterterm for the 1-loop power spectrum (4.12): 


rv rv y 

/ dr]i cir/2 Cafe(??,7?i) ^4'2(a3,??2)-VTb(ic,r/i). 
Jt] 0 Jr]l ^ 


(4.32) 


That this cancellation of IR-singularities must happen is almost obvious since any diagram that 
has a potential singularity (i.e. any diagram where the soft line attached to the loop is a v^di) can 
be identified with a contribution to (4.32). An explicit check is provided in appendix A.l. 


4.5.3 IPI diagrams: new counterterms 

After the subtraction of shift terms from IPI diagrams, what remains must be renormalized with 
the genuinely new higher order counterterms. These are made of locally observable operators 
integrated over the history of the short modes but evaluated at x (their displacement to x^ would 
correspond to higher order counterterms). At this level one must add a second order term to the 
equation of motion, corresponding to a cubic effective vertex. Prom (4.20), there are 18 possibilities 
for quadratic terms in c*-^. Note that even if ai = 02 the two soft legs of the amputated diagram 
can be embedded in different ways in a connected diagram and be integrated against different time- 
kernels, therefore no symmetrization is needed. When taking derivatives to construct counterterms 
did, the number of operators proliferates but some are degenerate. It follows that, for fixed 01^2 
and 771^2) possible c* vectors are 


diV‘^(j)aiV‘^4>a2, didj(j)aiV‘^dj4>a2, {«! 02, r/i O r/ 2 }. 


didjdk4>aidjdk4>a2 + {oi ^ 02 ,r/i o r/ 2 } 


(4.33) 


Thus there is a total of 15 independent, temporally nonlocal quadratic counterterms. It is shown 
in appendix A.2 that they are sufficient to fully renormalize the IPI diagrams. In fact, as expected 
from the structure of the leading vertices (4.4) and the requirement of momentum conservation (see 
the end of section 4.4), only a subset of these counterterms are allowed by the symmetries. The 
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correct prescription will be discussed in section 6. Also, because SPT does not generate vorticity 
(see subsection 4.3), an even smaller subset is needed to cancel UV-dependencies. 

In summary, we have demonstrated that the effective theory of large scale structure can be 
systematically renormalized, at the expense of introducing a plethora of temporally nonlocal coun¬ 
terterms. After discussing the implications of locality and momentum conservation in the next 
section, in section 6 we explore a more practical approach. 


5 Momentum conservation, locality and double softness 

As alluded to in section 4.3 momentum conservation and locality of the short scale dynamics ensure 
that the time evolution of short wavelength fluctuations can only lead to the generation of a longer 
wavelength perturbation suppressed at least by [24]. In the following we will first review the 
heuristic argument of Peebles, and next show that double softness follows from the equations of 
motion of the fluid system. 


5.1 A general argument for double softness 

First consider an almost homogeneous initial density po with small short wavelength fluctuations 
of scale 1/g, but no long wavelength fluctuation {S{k,r]o) = 0). Due to gravitational instability, 
these initial fluctuations collapse and form a clumpy distribution of matter p{x). We ask what is 
the typical size of a long wavelength fluctuation 6 {k, rf) after this process? Approximating the final 
distribution by a set of point particles of mass rUn at location x^, we have 

= (5.1) 

Pn 

where p is the mean final density. If the short-scale dynamics was turned off, the local mass 
distribution would be uniform; we denote it by a tilde: p{x) = p. There is a nonzero 5{k,r]) 
because this would-be uniform distribution has clustered, in this hypothetical uniform universe, 
to keep track of the matter that ended up in each clump, we divide the space into (possibly 
overlapping) regions with density pn{x). All of the mass in this region falls into the clump: 


mn= dP’y pn{xn + y). (5.2) 

Rn 

If there is no overlap Pn{x) is the same as p{x) (and hence uniform) within i?„. But if two regions 
Rni and Rn 2 overlap, then in the overlapping region Pmix) + Pn 2 {x) = p (and similarly for more 
overlaps). The typical size of each region is the extent to which mass elements are displaced from 
their initial position, namely, of order 1/q. Now we can write 


0 = -_ [ d^x-pe^^-^ = --y, f d%pn{xn + 

PJ PnRn 

(5.3) 

= 3 + ijedi - + 0{k^/q% 

P n 
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where d\^ and Qn are, respectively, the first and second moments of the mass distribution in as 
measured from £C„; 


d\ = 


' Rn 


d^y y* p{xn + y), 


Qn = f d^y + y)- 

Rn 


(5.4) 


The expansion of the exponential is justified because of the smallness of the size of the region 
compared to \/k. Comparison with (5.1) implies 


1 Jc^ 1^3 

6 ik,y) = ^ ^ 0 {k^/q^). 


(5.5) 


The second and higher order terms on the r.h.s. satisfy our general expectation about the depen¬ 
dence on k. As for the first term, the sum can be related to the center of mass position of each 
region: 

YdW’^-^- =Y^n{x^^ - (5.6) 

n n 

The difference (a;^'^ — Xn) can only be caused by momentum transfer among different regions. If 
initially there existed only modes of wavelength ~ l/g, there is no coherent momentum transfer 
over much longer distances, and the sum over any patch of size much larger 

than 1/q is negligible. Hence, the sum in (5.6) is nonzero to the extent that the exponential varies 
over any such patch, and therefore it must be suppressed at least by one factor of k. This shows 
that all terms in (5.5) are suppressed by at least two powers of /c.® 

In the absence of long wavelength initial perturbations the location of the clumps Xn is effectively 
random at large scales, and the sums in (5.5) become identical to random-walks. Correlating two 
such contributions and averaging over the short modes results in 

.^stochastic (^ ^ 0) OC A: . (5-7) 


This is the generic asymptotic behavior of stochastic contributions to the long-wavelength power. 

In the presence of initial long wavelength perturbations the above argument can be modified 
as follows. Now the clumpy universe should not be compared with a homogeneous universe, but 
with one in which the long modes are kept but the short modes are absent. Hence, p{x) is smooth 
but nonuniform. All other steps remain the same, though one should bear in mind that x^ and 
both depend on and can get largely displaced by an initial long mode 5i(fci); the leading 
dependence which goes as 6 i{ki)/ki cancels in the difference in (5.6). Nevertheless because of the 
remaining dependence of both x^ and the shape of each region i?„ on the initial long modes, (5.5) 
has nonzero correlation with the initial (5i(fci). The non-stochastic diagram (3.1) (equivalently 
(3.8)) is one perturbative analog of such contribution. 

Finally, a similar argument implies that the momentum density 

7r(fc, p) = (5.8) 

^ n 


starts from order k. 

®Our argument slightly differs from that of [24] in that the first term on the r.h.s. of (5.5) is absent there. However, 
this seems to be a nontrivial part of the argument. 
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5.2 Double softness in pertnrbation theory 

The double softness property is not manifest in the formulation (2.1) in terms of 5 and v. However, 
it becomes manifest if we eliminate v in favor of the momentum density tt = (1 + 5)v, in terms of 
which the continuity equation becomes linear (4.17) and the Euler equation becomes 

re'-+ 'H'k^ =-{1 +5)di(t)-dj . (5.9) 

Writing 

ddicj) = dj{dj(j)di(j)) - ^di{dj(j)dj(j)) (5.10) 

shows that the non-linear terms in this equation all have an overall derivative. Moreover, taking the 
divergence of this equation and using dire' = —6 to obtain an equation for 6 makes all interaction 
vertices explicitly second derivative. In appendix B it is shown that these equations evolve ingoing 
hard lines into a single soft answer for 7r(fcout) and a double soft answer for 5(fcout) when kout ^ 
{/cj}in. The same conclusion is shown to hold when there are ingoing soft lines with momenta of 
the order of kont in addition to the hard lines. It then becomes clear that the softness property is 
preserved if counterterms of the form djT'^ are added to equation (5.9), where is made of local 
observables. This ensures locality and momentum conservation. 

6 A practical remedy 

So far we learned that a systematic renormalization of the two field (<5, 0) formulation of cosmic 
fluid is possible if one uses nonlocal in time counterterms and add corrections to both continuity 
and Euler equations. However, there are three reasons to seek an alternative approach: 

I One may not be interested in the velocity statistics, and even if she is, there are various ways 
to define velocity field and the one used in section 4 might not agree with the velocity field 
obtained from an actual survey or simulation. 

H More importantly, the momentum conservation constraint is not straightforwardly imple¬ 
mented. 

HI The nonlocal in time counterterms introduce additional complications, since instead of a single 
time-dependent coefficient one needs to consider a multi-variable function of insertion times 
for every new counterterm. 

We first show in section 6.1 that one can use the momentum formulation to eliminate velocity. In 
this language, only the Euler equation receives corrections and not the continuity equation. Also, 
now there exists a straightforward prescription to generate all counterterms and only the ones 
allowed by the symmetries, including the conservation of total momentum. The observed velocity 
on the other hand must be treated similar to biased tracers; it can differ from the one calculated 
in any perturbative scheme by a set of locally observable counterterms whose coefficients must be 
empirically determined. This resolves the first two problems. Then in section 6.2 we show that all 
nonlocal in time counterterms can be reduced perturbatively to a set of local operators, addressing 
problem HI. This allows us to give a complete and more manageable list of counterterms. 
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6.1 Formulation in terms of the mass-weighted velocity 


If we use the density 6 and momentum tt as the primary variables, the continuity equation reduces 
to the linear equation (4.17), and does not require any counterterm (it is not modified by smoothing) 
[2, 18]. If the correlators of 6 (and therefore of 6 ) are renormalized, then so are the correlators 
of V • TT and vice versa. The same can be said for (j) and 6 , given that the Poisson equation is 
also linear. The problem then reduces to solving the only non-linear equation, namely the vector 
Euler equation (5.9). If one is interested in the correlators for 6 , it is convenient to decompose the 
momentum into a divergence (// = V • tt) and a vorticity (iz = V x tt), as discussed in appendix B 
(see (B.l), (B.2) and (B.3), see also [31]) . The Fourier space version of these equations is collected 
in appendix D. Notice that all nonlinear terms on the right hand side of the Euler equation (5.9) 
are total derivatives of some tensor that is not proportional to Sij and hence all contribute to the 
curl equation. In particular, all non-trivial solutions have i ',6 ^ 0. Using the continuity equation 
(4.17), one can get rid of jj, and solve the remaining equations for 6 and iz. When using momentum, 
one can therefore compute all correlators of 5 without including any counterterm in the continuity 
equation. We will see now that this is still the case when using a specific definition of the velocity. 

In Standard Perturbation Theory (SPT), where Tij = 0, there exist a convenient change of 
variables for which the curl of the vector Euler equation is trivially satisfied even for non-trivial 
solutions (5 7 ^ 0. From (5.9), it is clear that a necessary condition for this to happen is that the curl 
is not sourced by the term (1 -|- d)di4>. Upon dividing the whole equation (5.9) by (1-1-5), this term 
reduces to the gradient of a scalar which disappears after taking the curl. This suggest to use the 
variable v.,^ = tz/{\ + 5)^ which, as it is well-known, indeed decouples the vorticity V x from 5 
in SPT. We will refer to v-,^ as “mass-weighted” velocity. 

More explicitly, let us define the divergence and curl of the velocity (we raise and lower all 
indices with 8 ij) 

e^ = V ■v^ = divi, wi = V = (6.1) 

from which it follows diw]^ = 0. The equations of motion then take the schematic form (see 
appendix D for the explicit form) 


8^-5 -\- 

dr9^ + ne^ + 

driv-K + T-Lwir 


—adj^d + dw ■ Wn8 , 







w. 




w^wt 


+ d X 



( 6 . 2 ) 

(6.3) 

(6.4) 


We can now perturbatively solve for 0-,^ and Wt^ in terms of 6 using (6.2) and (6.4). This yield, 
order by order, a single differential equation for 6 . Let us see explicitly how this works to quadratic 
order. The vorticity is sourced by (9 x [dr /(1 -|- 5)] and therefore starts quadratic in 5 perturbations. 
We can hence neglect it in the following formulae, keeping in mind that it will contribute to 5 at 
the next order, 0(5f). We now solve the continuity equation 


=-5 + + o{5^), 


(6.5) 


and substitute it into the divergence of the Euler equation to obtain a nonlinear equation purely 
in terms of 5: 


s + m- = {dr + n)di{6^6) + ^^{^6^6) + o{6^). 


( 6 . 6 ) 


24 





The linear part of this equation is a second order differential equation which gives a retarded 
Green’s function, using which one can perturbatively solve for the evolution of 6{t,x) in terms of 
the initial condition 5i(»), without ever mentioning 0.,^. Obviously the resulting equation for 6 now 
has interactions of all orders in 6 . On the other hand, the original system (2.1) contains only cubic 
interactions and quadratic mixings between intermediate 6 and 6 lines. What the above procedure 
of solving for in terms of 5 means is that all intermediate 6 lines are contracted to points, thereby 
generating arbitrarily high order interactions of 6 . 

This formulation justifies a common approach in the EFT literature which is to use 5 and 9 
as variables but to consider counterterms only in the Euler equation while keeping the continuity 
equation intact. It is equivalent to renormalizing the equation (6.6), as any conventional nonlinear 
theory of a single field <5, and to redefine the composite operator 9^^ to absorb the counterterm V ■ jc 
of the continuity equation, so that the redefined 9-,^ is given by the same relation (6.5) at all scales. 
However, doing so can in principle introduce vorticity because V • jc is not necessarily Laplacian 
of a locally observable scalar. Therefore, the counterterms of the (vectorial) Euler equation needed 
to renormalize the SPT loops are no longer pure gradients. This was the case in [21], where 
the quadratic counterterms dj{didk 4 >djdk 4 >), which is not the pure gradient of a locally observable 
scalar, was needed to renormalize one-loop bispectrum of 6 . 

In summary, one can use the standard perturbation theory approach to calculate various corre¬ 
lation functions of 6 . The counterterms that are allowed to be added to the vector Euler equation 
are now fixed by momentum conservation to be of the form 

{l + 6 )-^djT^^ (6.7) 

where in order to do systematic renormalization must be allowed to depend on observable 
quantities along the flow, and hence be nonlocal in time. 

As for the statistics of the observed velocity field, the difference with that of perturbation theory 
(say Utt) can be parameterized by a set of counterterms 

Wbs ^TT T Wt) (6'^) 

where Vct is a locally observable vector. It has a non-stochastic piece which by symmetry must start 
from first order in derivative counting, and a stochastic piece which according to the arguments of 
appendix B.3 is generically 0{k^), leading to 0{k‘^) contribution to the vorticity power spectrum. 
Hence, the treatment of velocity is similar to the biased tracers. In particular the most general Vct 
can be expressed locally in time [8]. 

6.2 Local formulation 

It was argued in the last section that a generic counterterm can be constructed out of a tensor . 
The tensor is made of observable quantities integrated along the flow: 

T*^(a;,r)= [ {fir„}ifa-a„-(r; {r„}) • • • na’;^’*(sfl(rn), r^) • • • , (6.9) 

Jo 

where we switched to the conformal time r. In perturbation theory, any observable can be solved 
in powers of the initial field 6 i{x), in particular, 

n(a;,r) = {x) + D‘^{t)U^^\x) + ■ ■ ■ (6.10) 

where D{t) is the growth factor and we are using the accurate approximation of replacing a(r) —>■ 
D{t) in going from matter dominance to ACDM. This expansion contains displacement terms to 
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take into account the bulk flow induced by long wavelength initial perturbations. As such, some 
of the second and higher order terms in the sum become singular in the presence of a finite di{ki) 
with vanishing momentum ki —)• 0. These terms are responsible for shifting the argument of the 
fields on the right hand side (r.h.s.) to their initial (Lagrangian) fluid position 

z = a;fl[0; £c, r]. (6-11) 

After shifting the arguments a; —>■ z there remains no infra-red singularity on the r.h.s. as expected 
on symmetry grounds. We denote the operators thus obtained Lagrangian operators. In terms of 
them [21, 8] 

Uix,T) = i?(r)ng(z) + D\t)u\11{z) + ... (6.12) 

Substituting this expansion in (6.9), makes it possible to perform all time integrations to obtain 
an expansion of various products of Lagrangian quantities with coefficients that generically depend 
on (and only on) the final time r. The whole expression is acted on by two Eulerian derivatives: 


A 

dx^ 


_^ 

1 -b (5 dx^ 


{rrin} 


(m„) 

Igr 


(^) 


(6.13) 


It was shown in [8] that there is a one-to-one correspondence between the Lagrangian quantities 
and the locally measurable (Eulerian) operators. Hence, we can re-express (6.13) in terms of those 
Eulerian operators. A basis for the latter can be obtained by combining various convective time 
derivatives (4.10) of didjcj). Starting from 

fl\f{x,T)=didj4){x,T), (6.14) 


the convective derivatives can be combined in such a way that a set of operators with increasing 
order in perturbation theory is obtained. Suppressing the indices and the argument {x, r) these 
higher order operators can be computed recursively: 


n(»*) 


(n(ri-i))' _ (n- i)n(”-i) 


(6.15) 


where prime denotes convective derivative with respect to logD (= 77 in a matter dominated 
universe). At leading order these operators match those that appear in the perturbative expansion 
of didjcf)-. 

nW(a;,r) = 14"(r)ng(z) + • • • (6.16) 

At lowest derivative level an Eulerian basis for two-index tensors can be constructed as follows. 
Eirst, a two index tensor is made out of a (matrix) product of a (possibly empty) sequence of Lt^’^^’s. 
Eor instance, djj, are two examples. This tensor can be multiplied by an arbitrary 

scalar constructed out of traces of products of ^ij \ like The trace Trp^”^] with n > 1 

must be excluded since it is expressible in terms of lower order operators. The first few orders are 
(in matrix notation): 

l^t lTr[n(^)],n(^) (6.17) 

2“'^ lTr[(li(^))2], l(Tr[li(^)])^ H^^^Trp^^)], (n^^))^ H^^^ 

3"^ lTr[(n(^))^], lTr[(n(^))2]Tr[n(^)], l(Tr[^(^)])^ ITrp^^^H^^)], •.. 


Note that the matrix must be symmetrized. To obtain a list of counterterms, one substitutes 
any such tensor into the Tij of (6.7). Let us close with a few clarifications. 
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• This is a complete basis at lowest derivative level. In particular, all locally observable opera¬ 

tors made of the velocity field are also included. For instance, in terms of a rescaled velocity 
potential we have 

1l[f = didm4>djdm<P + ^{didjCp - didjcpy) + 0{6^). (6.18) 

Clearly, at higher orders didjcj) and didjcj)^ will be insufficient for a local description of coun¬ 
terterms and higher order convective time derivatives will be necessary. However, given that 
at order we need up to the first example appears at the quartic level. 

• The tensors ri^”^(a;,r), and hence the counterterms in (6.7), are easily calculable in momen¬ 
tum space by combining the F and G kernels of SPT. For instance, 

n!f(fc,T)= (6.19) 

' ' n=l 


where 6n{k) is the order SPT solution, and 




n=l 




(n - l)6n{k) 


and so on. 



{k — p) ■ p p^p^ 

\k — Pp IpP 


Sm{k p^dji—m 



( 6 . 20 ) 


• One should bear in mind that although didj4>{x, r) is a total derivative, neither its convective 
time derivatives nor the 2 :-space matrices Il\j\z) are necessarily total spatial derivatives of 
any kind. 

• Although, the basis (6.17) for tensors is non-redundant, taking the derivatives leads to degen¬ 
eracies. For instance, at second order there are 5 tensors but only 3 indepedent counterterms. 

• Any tensor leads to an infinite series of counterterms when the factor (l-l-d)”^ is expanded 
in (6.7). The relative coefficients of these counterterms are fixed so as to insure double softness 
of 6 when the ingoing fields are hard (a fact that is verified in appendix B using the momentum 
equations). 


7 Conclusions 

In this paper, we have introduced a systematic renormalization procedure for the perturbative 
description of Large Scale Structures, which is analogous to the one commonly used in quantum 
field theory. One main advantage is that new counterterms can be straightforwardly read off 
amputated time-evolution diagrams, without the complicated mixing with lower-order counterterms 
that plagues the standard treatment. In this approach, correlators of both 6 and 6 are renormalized 
at each order. We have also provided for the first time an explicit prescription for the allowed 
counterterms to all orders in perturbation theory. In passing, we have justified the commonly used 
mass-weighted volcity formulation in which counterterms appear exclusively in the Euler equation.® 

®This had been checked so far explicitly only for the one-loop power spectrum [20] and at linear order in the 
leading counterterm to the renormalized velocity [18]. 
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We have also presented a rigorous derivation of the double-softness of 5 and V ■ tt in perturbation 
theory, namely the well-known fact that 6{k) ~ V7r(fc) ~ k‘^ on large scales, A; —>■ 0. 
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A Second order counter terms at 1-loop order 

This section is devoted to a detailed calculation of second order counterterms in the leading loop 
order. This is sufficient for a 1-loop calculation of the bispectrum. However, the resulting coun¬ 
terterms renormalize all reducible diagrams where the vertex corrections of Fig. la and lb appear 
as a sub-diagram, which is the goal of systematic renormalization. For this the counterterms are 
unavoidably nonlocal in time since the ingoing lines into this sub-diagram can be of different orders 
and hence have different time-dependence. This is the practical manifestation of the “long memory 
effect” discussed in section 4.1. 

Using Feynman rules presented in section 2, the 1-loop correction to the vertex function can be 
read off diagrams of Fig. la and Fig. lb as 


Cla°a2 (h, hi, ^2 i ^, fcl, ^2 ) + C'a'aj^aa (h, hi, h2 i ^, ^1 ,^2 ) 


/ ki ^ k2 
+ hi ^ h2 
\ ai 02 


(A.l) 


for which the diagram.1 and diagram.2 can be calculated through performing hard momentum (loop 
momentum) integrations of the following equations. At this step time integrals are left untouched 

Caa^a2iv,VuV2-,k,ki,k2) = 4 X gfm{v2,0)ghi{vi,0)gbe{v,vi)9cd{v,m) 

/ la,b,c{k,ki +q,k 2 -q)ldj,a 2 (^2 - q,-q,k2)7eaih{ki +q,ki,q) 
Jq 

(A.2) 


Caai\ 2 ('n,Vi,V 2 -,k,ki,k 2 ) = 8 X gbi{v,0)ge,m{Vi,0)9cd{v,V2)gfh{V2,Vi) 


/ labc {k, -q, k + q) Jdfa 2 (k + q,ki+ q, ^ 2 ) Ihea^ (ki + q, q, ki) 

Jq 


(A.3) 

and summation on repeated indices is presumed. It should be noticed that the counterpart of 
diagram. 1 which can be found by exchanging ki ^ k 2 ,rji h 2 )fli 0 - 2 , would just multiply 
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ki <-> k2 

m <-> m 

ai o 02 



fci O fc2 \ 

I 

Oi 02 / 


(a) Amputated Diagram. 1 


(b) Amputated Diagram.2 


Figure 1: 1-loop corrections to the cubic vertex function 


the contribution of diagram.1 in Ca~^^°a 2 by a factor 2. However, in the second diagram, (1 2) 

conversion produces a new momentum and time dependence. 

The results of the above integrals are 


= cr’^ 0 {r] - r]i) e{r]-r] 2 ) 


IplD „2r, , 1t.(2) (l/2)(-,7+5r,i) , 

^ a,ai,a2 ^ ' a,ai,a2 ^ ' 


IpW ^ g(l/2) (-»?+5 t; 2) _|_ IpD) ^ g(5/2) (tu+?72)-3t; I (A. 4) 


C'a,a;^a2(^;^D^2;??,??l,??2) = 0{V - m) 0{r]2-Vl) 

12t.(1) 2r, , 2t.(2) (1/2) (-,7+5t;i) , 

^ ^ a,ai,a 2 ^ ' a,ai,a 2 ^ ' 

2p(3) g(l/2)(-7?+5»,2) , 2t.(4) (5/2)(r,i-r,2)-t2,? 1 /ac) 

^ a,(11,0.2 ^ J- a,0.1,0.2 ^ f [^.OJ 

in which and are originating from 1-loop vertex diagram 1 and diagram 2 respectively. 

A.l Shift terms 

As discussed in section 4.1 promoting the time integration to be along the flow identifies several 
higher order counterterms as the displacement a lower order one. These shift terms can be isolated 
by looking at the infra-red limit when one of the ingoing momenta becomes much smaller than the 
other. In other words, the two vertices v^diO and v^diS with soft velocities lead to IR-divergences 
as 1/k in the momentum of 'I '2 ~ 0 = V.v. In the case of 1-loop second-order counterterms, the 
first order shift of a first order counterterm (4.13) automatically renormalize these “IR-divergent” 
cutoff dependence. It is shown in section 4.5.2 that these non-original counterterms are 

p n V 

/ drji dT]2 Cab{ri,rii) (A.6) 

Jrjo Jrji ^ 
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Now using the 1-loop correction of the propagator 
Cab{r],Vi) = -(^v 


f 

3 

-4 

5t7 

' ^ A ' 

5t7i I 

) 

1 - 

25 

-8 

5 

e 2 + 

5 25 

9 3 

5 5 

e 2 


the above correction to the equation of motion can be written in Fourier space as 
fv rv 

/ drji drj2 Caaia2iv,m,'n2-,k,ki,k2) '^2{k2,'n2)- 

Jt]o J Vo 


(A.7) 


in which 

Caaia2iv,m,'n2',k,ki,k2) = -Cr'^ 6a^2 


2 k'^ {ki.k2) 


kl 


rv rv 

/ drji / drj 2 9{r] - rji) e{rji - 7 ) 2 ) 
Jrjo Jrjo 



1_ 

1 tol 1 


1 

L_ 

5771-77 

e 2 

1 

5 

5 

d di 

5 5 

aai / 


(A. 8 ) 

Now let us examine the IR-limit of the corrections to the vertex function. Using the following 
identity between step functions 


9 {r] - r]i) 9 {r] - 7^2) = 9 {r] - r/i) 9 {r]i - 772) 9 {r] - 772) 0(772 - 771) 

one readily finds that IR-limit of diagram. 1 contribution to the vertex correction as 

{ki.k 2 ) n ^ n 


(A.9) 


+ {I ^ 2) = - 2 x a 5a^2 


^2—>-0 


rv rv 
/ drji / dr]2 

J T}' Jr}' 

9 {r] - r]i) 0(771 - 772) 9 {r] - 772) 0(772 - 771) 




_ 6 _ 

25 17| 




a ai 


-6 -12 
25 175 

0 4 


5V1-V 

e 2 


a ai 


(A.IO) 


where the symmetry factor 2 comes from (1 -H- 2 ) or more explicitly (fci -f-)- k 2 ,r]i o 772,01 -H- 02 ). 
And 


lim„'5C'aaia2(^>^l’^2;7/,7/l,02) "b (1 ^ 2) 


— 6 , 


fc2—>-0 


0(77 - 772) 0(772 - 771) 
-0(77 - 771) 0(771 - 772) 


a2 2 


-12 -24 

25 175 

0 — 

^ 25 


_3_ 

5 25 


fcn 


o2»? , 


J a ai 


rv 

rv 

/ 

d77i 

/ dr]2 



Jv' 


12 

24 


+ 

25 

0 





2^ 

a 




J d di 


-28 52 

5 25 


Ovi-v 

e 2 


e 2 


J a ai 


(A.ll) 


Again the third line is the soft limit amplitude of the inverted version of Diagram.2. Adding up 
these contributions together one recover the IR-divergent vertex correction already anticipated by 
promoting time integrals in propagator correction integral, , to an integral along aifl, which 
leads to (A.8). 
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A.2 The new counterterms 


As discussed in 4.5.3, after the subtraction of shift terms from IPI diagrams, what remains must be 
renormalized with the new higher order counterterms. For every second order term (counter-term) 
in the equation of motion there is an associated cubic vertex in Feynman diagrams so we use 1-loop 
second order counterterms and 1-loop cubic vertex correction interchangeably. Using Eq. (6.2), all 
second order counterterms are of the following form 



dm F\k,ki,k2) 4'ai(^i,m)^a2(^2,h2)- 


The general structure of E* kernels can be found using Eq. (4.20) as 


F^{k,ki,k2) = k.ki 


F^{k,ki,k2) 

F\kMM) 


(fc.fci)(fci.fc2) 

kl 

{k.k^){ki.k2f 

U2 u2 
^2 


F\kMM) 

F\k,kiM) 


F\kMM) 


{k.k2) 

{k.k2){kl.k2) 

kl 

{k.k2){ki.k2f‘ 

U2 jL2 


(A.12) 


(A.13) 


As already emphasized, although the basis 4.20 for tensors is non-redundant, taking the derivatives 
leads to degeneracies. It can be shown that the above structures are linearly dependent: F^ — F^ = 
F^ — F^. Hence, we eliminate E® — E® and define 


E® 


1(F' + F«). 


(A.14) 


Having found all possible momentum dependencies, the next step would be listing the time kernels. 
There is a unique basis for all possible time kernels of 2-nd order counterterms in the equation of 
motion. We denote them by k, given by 


kl = ka = ^ g(57?2-r,)/2 

k4 = e5()?i-’?2)/2-|-2»; _ g5(??i-»;2)/2+2r; _ g5(?7i+772)/2-3r; 


(A.15) 


Every time kernel has three indices associated with indices of the one outgoing and two ingoing 
legs respectively. As every index can be either 1 or 2 there are 2^ = 8 distinct counter terms which 
are listed below 


'y ^ Ari2iE* — ~l~'l25^^ ^ — 6(4ki — dka — 7k5 -|- 7kg)0i2 — (24ki — 4k3 -|- 63k4 -|- 42k6)02i] 

i 

+ [3(44ki -I- 25k2 - 74k3 -|- 3k4 -|- 2k6)02i + (132ki - 57k2 - 200k3 - Okg -|- 6k6)0i2] 

- [6(13ki - 13k2 - 4k5 -|- 4k6)0i2 + (78ki - 13k3 -|- 36k4 -|- 24k6)02i] (A.16) 

^ ^ K222^^ ~ [(“74^1 + 25^2 + 44^3 + 2/^4 + 3^6)021 + 

i 

(—74^1 -|- 44^2 “1“ 25^3 “h 2^5 3^6)012] 

— ——F^ |^(—13^1 + 78^2 + 24^5 + 36^6)012 + (—13^i + 78^3 + 24^4 + 36^6)02i] 

i.ZtO 

(A.17) 

V [02i(129ki + 70k2 - 199k3 + 33k4 - 33k6) + 30i2(43ki - 43k2 + llkg - llkg)] 

875 
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Figure 2: Hard ingoing lines combining into a soft line. 


2 

+ [012(129ki — 199k2 + 70k3 + 33% — 33%) + 3021 (43ki — 43% + 11% — 11%)] 

+ + F"^ + F^) [0i2(-611/ti + 611% + 48% - 48%) 

o75 

+ 02i(-611% + 611% + 48% - 48%)] (A.18) 

y~] Kl^iF^ = -^Fi [0i2(6ki - 16% + 15% - 3% - 2%) + 02i(6ki - % - 3% - 2%)] 

i 

+ ^F2 [021 (6ki + 15% - 16% - 3% - 2%) + 012(6^1 - % - 3% - 2%)] (A.19) 

.Zo 

= 7^(7^^ + F"^) [0i2(16ki - 6% - 15% + 2% + 3%) + 02i(16ki - 16% - 3% + 3%)] 

I 

(A.20) 

7^211-^* ~ [012(~2ki + 2k 2 + ^5 ~ ^e) + 02 i(“2ki + 2k 3 + % — Ke)] (A.21) 

i 

in which 0i2 = 0{r] — rii)6{r]i — 772 ) and 02i = 0{ri — r]2)0{r]2 — ill)- Moreover and 

Si 7^212-^* could be simply read off from K\ 2 iF'^ and K 22 iF'^, respectively, by the following 
replacements 

K2 %) % ^ ^5) F^ F'^ , F^ F'^, 012 ^ 021 (A.22) 

Note that the counterterms for the Euler equation, which are the ones with the first index 2, 
are all proportional to 

F, + F 2 = k\ or = (A.23) 

Ki ^2 

which in real space translate into and V^(clic)j%)^, respectively. Hence, as expected they are 
Laplacians of locally observable scalar quantities and the corresponding force in the Euler equation 
is curl-free. 

B Softness properties of perturbation theory 

In this appendix we give a perturbative and iterative proof of the single softness of momentum 
7r{k) and double softness of density contrast S{k) which result from evolution of short wavelength 
initial modes - also known as decoupling. The proof is broken into several steps. 
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B.l Only hard ingoing modes 


Let us first list the equations and constraints to be satisfied by <5, = V x tt, and ^ = V • tt, 

namely (4.17) and the divergence and the curl of (5.9) 

S = -V-7v, (B.l) 

6 + T-Ld — = didj{di(l)dj(j)) — ^\/‘^{dj(l)dj(j)) + didj{v''7r^), (B-2) 

= -e^^’^dj[dm{dkcl)dm(p) + d^iv’^TT^)] (B.3) 

where we used 

TT 

TT5 

and 

VV = ^'H‘^5 (B.5) 


to simplify the last two equations (and the subscript on is dropped to avoid clutter). Suppose 
we start from a set of initial short wavelength (hard) modes. These modes evolve and combine 
according to the above equations to generate second and higher order perturbations. The softness 
properties follow from; 


1 . At any level when, as in figure 2a, hard modes combine into a soft 5 through the vertices of 
(B.2) 

Hv) = [ dr]'Gs{v,v')[didj{di(l)sdj(l)s) - ^V'^{dj(l)sdj4>s) + (B.6) 

dvo ^ 

the result is automatically of second order in the outgoing momentum, i.e. kf. (Here and in 
what follows we use the more compact subscripts I and s to denote long wavelength (soft) 
and short wavelength (hard), respectively.) 

2. The same applies to V • vr; because of the constraint (B.l), and (B.3) ensures in a similar way 
that V X TT; = 0{kf). It follows that tt; = 0{ki). 

3. If the hard modes combine into two long modes at a lower order which then evolve into one 
long 5i (as in figure 2b), we have 

5i{'q) = f dr]'Gs{r],r]')[didj{di4>idj(j)i) -^V‘^{dj(t>idj(l)i) + didj{vl7rj)]rj:. (B.7) 
dm ^ 

The softness of the final bi can be established by induction as follows. 

4. If all long modes on the r.h.s. are the immediate result of combining short modes, we have 
from 1 and 2, and (B.4) and (B.5) 

di4>i = 0{ki), vi = 0{k^), ni = 0{ki), (B.8) 

which ensures that 6i on the left of (B.7) is at least 0{kf )d^ This in turn ensures via (B.l) 
that V ■ TT/ = 0{kf) at this level, and one can similarly verify i' = 0{kf). 

5. Therefore a stronger softness property than (B.8) holds when the soft modes on the r.h.s. of 
(B.7) result from evolution of two soft modes that themselves evolved from hard modes, and 
the iteration can be continued ad infinitum. 


^'^Since (B.4) is nonlinear and has no overall derivatives a priori one would expect vi = 0{kf). However, we will 
show in appendix B.3 that at lowest derivative order the (standard) perturbation theory results in vi = 0{ki). 
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B.2 Soft and hard ingoing modes 

In the presence of both soft and hard initial modes, double softness can be established by dividing 
the quadratic vertices of the perturbation theory into different categories. First, every vertex that 
results in a soft outgoing field either combines two short modes or two long modes, e.g. 

[d(t)d4)\i = dipsdips + d4>id4>i. (B.9) 

In the presence of initial soft modes we need to make a further distinction of whether a long mode 
is purely made of those initial long modes (indicated as d(pi-i), or if it contains initial hard modes in 
addition to a (possibly empty) set of initial long modes (indicated as dcjii-si)- Thus, the first vertex 
in (B.2) is decomposed schematically as 


dd[d(l)sd(t)s + d(j)i.id(l)i.si + d(j)i-sid(l)i-,si]- (B.IO) 

Note that d(j)i-id(pi-i is excluded since we are only interested in the case where there are some hard 
initial modes. Figure 3 shows one example for any of the above three terms. 

1. The first term in the above equation is similar what we considered in point 1 of the previous 
section, where the two external derivatives ensured double softness. However, there is an 
important difference: the short modes can now contain initial long modes (see figure 3a). 
One may therefore worry that d(j)s may contain negative powers of ki due to the displacement 
of hard modes by the soft modes. That is, terms in perturbative evolution of the form vi-i • 
where vi-i = 0{kf^). These IR singularities indeed exist in any short wavelength quantity 
that evolves from initial hard and soft modes. They correspond to the motion of short 
modes with the bulk flow of the long modes and their presence is ensured by the equivalence 
principle. However, as shown more explicitly in appendix C whenever a collection of hard 
modes combine into a soft mode, adding any number of initial soft modes to such diagram 
does not lead to any negative power of ki. The singular displacement of hard lines with 
momenta {qi} by a soft line of momentum ki that are proportional to ki ■ Qi/kf add up to 

~ 1 = 0{k^) at any moment. Hence, the first term in (B.IO) leads to a double 
soft outgoing 5i (and as before a single soft tt;). 

2. The second and third terms can also be shown to lead to at least double soft results using 
induction. First suppose there is no initial long modes in dcpi-si- Then by the argument of 
the last section d(j)i-s = 0{ki). And d(j)i-i = 0{k'f^) from (B.5). So the resulting outgoing 

= 0{kf) because of the overall derivatives. One can then continue the induction provided 
that the same argument applies to the dd{v7T) vertex of (B.2). 
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3. Decomposing vtt into short and long wavelength contributions we get 

[vTt]i = VsTTs + Vl-sm-l + Vl-slT^l-sl + Vl-lTTl-sl- (B-H) 

Assuming that up to some order in perturbations 

n-si = 0{h), 5i,si = 0{kf), vi.si = 0{kf), (B.12) 

and using the fact that vi-i,Tri-i = 0{kf^) and 6i-i = 0{kf), we see that the last two terms 
on the r.h.s. of (B.ll) preserve these properties (when considering the vertex dd{v7r)). The 
second term clearly violates them and the first term does so too: expanding 


VsTTs =Vs[SsVs + (1 + S)iVs + SsVi] 


=Vs[SsVs + (1 + S)iVs + SsVi-si + dsVi-i], 


(B.13) 


the last term is seen to be 0{ki ^), leading to 0{ki) contribution to 5i from the vertex dd{vsT^s)- 
However, the sum of the two dangerous terms preserves the condition (B.12). Substituting 


■Ki-i = {I + 5)i-ivi-i (B.14) 

in the second term on the r.h.s. of (B.ll) and using the fact that (1 + S)^si = ^i,si (since the 
subscript s indicates a nonzero set of initial short modes) we obtain for the sum of potentially 
dangerous terms 


^'Z;z[(l + S)l-lVl.sl + = Vl-l {[(1 + 6)v]l.sl “ (1 + ^)l-,slivi-,l + Vl-,sl)} 

(B.15) 

= vi-ilni-si - + vi-si)] = 0{ki). 

Therefore the sum of all vertices preserves the condition (B.12), and double softness follows 
by induction. 

Note finally that all of the above arguments generalize to the case where counterterms of the 
form djT^^, with a locally observable tensor which by definition can never have an IR singularity 
(corresponding to a negative scaling with ki), are added to the Euler equation (5.9). This seems to 
be the most general way of modifying the equations of motions that is still compatible with locality 
and momentum conservation. Intuitively, the short scale dynamics can generate an effective force 
on a volume element of the fluid only due to its momentum flux into that volume. 

B.3 Velocity field 

The final point to be discussed is the softness property of the mass weighted velocity field defined 
as u = 7r/(l + (5). Using mass conservation we get 

0 = V-r; = -j-V-((5u). (B.16) 

Using this equation in the Euler equation leads to 

5 + n6- = n{dr + n)dii6di<p,) + (b.it) 

where we neglected vorticity and set v = —T-LV(j)y. From these equations it follows that double¬ 
softness of a soft outgoing 5 evolved purely from initial hard modes implies double softness of 9 . 


35 


To show this one needs to argue that the second term in (B.16) which naively is only single soft is 

actually double-soft when the perturbative solution for 5 and 0 are substituted. This follows from 

the fact that the second vertex on the r.h.s. of (B.17) is automatically double soft and if 5 is to be 
double soft the first term on the r.h.s. must also be double soft by itself. 

If we use 8 and v as primary variables - as we did in our systematic renormalization - and 

add counter-terms proportionally to SPT cutoff dependence, the resulting renormalized velocity 
inherits both curl-freedom and single-softness from SPT kernels. This renormalized velocity of 
course differs from On the other hand, if we use and add counterterms as to the 

momentum equation (5.9), they modify (B.17) as 5j[(l -|- Expanding this leads to 

vertices which are not explicitly second derivative. For instance, taking = 5^^ {dmdn4‘)‘^ we get 
an associated cubic vertex 

- di[8 di{dmdncl)f]. (B.19) 

Taking all ingoing modes to be hard, this vertex gives rise to an outgoing 8 which is only single soft. 
In this formulation, the full result becomes double soft (as it must according to the arguments of 
sections B.l and B.2) because the above single soft term cancels with the one arising from evolving 
a second order solution coming from the leading part of the counterterm {dmdn4‘)‘^ to third 
order via the vertex di{8v^). Hence, the double softness of 8 does not imply that of 9-,^ anymore 
because the second term on the r.h.s. of (B.16) is generically only single soft after the addition of 
counterterms. Note that one can redefine velocity to absorb these single soft terms, but the new 
velocity is no longer given by 7 r/(l -|- 8). 


C Cancellation of IR singularities 


Consider the following time evolution diagram with ingoing momenta {fci, 91 , 92 } and outgoing 
momentum k, and such that /si ~ A; <C 91 , 2 - 


k 



(C.l) 


The potential IR singularity arises in this diagram from the vertex v ■ V'I'a, where v realizes the 
soft ki mode and 'ka a hard one, leading to ki ■ qi/kf. Focusing only on this vertex we have for 
the outgoing field 

rv 

{^a(A:,9)}lR = / dr]igabir],r]i) 

^2^^(fci,92)^i^^(9i,92)^d ^( 92 , 91 ) + 9i ^ 92 ]- 

^^One should also consider the case when the hard modes combine into two soft modes which later combine into a 
final 5. In this case too the second term on the r.h.s. of (B.17) 

(vi-sVi-s) (B.18) 

is double soft because at least vi-^s = 0{kf). 


/ dr]2[Ybcd{k,qi + fci, 92)gee(9B 92 ) ,2 

Jr?o % 


(C.2) 
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Now using 


'^d{q2,m) = 9df{m,V2) ^f{q2,v2), 


(C.3) 


and the fact that 'jbcdi^, qi + ki, 92 ) = 'Jbcdik, Qi, 92 ) + 0{ki) in the limit ki <C Qi, we obtain 

rn pi 

{'^a{k,q)}iR= dr]igab{v,Vi) / dr]27bcd{k,qi,q2)gce{vi,V2)9df{vi,m) 

^ VO VO 

(C.4) 

{92,92)- 

The expression in the square brackets is ki ■ {k — fei)//c^ which is 0 {k^), thus the potential IR 
singularities coming from the v ■ vertex cancel in this diagram when we sum over the two 

possible attachments of the soft line. The full diagram will therefore have the same softness property 
as the one without attaching the soft line. 

The above cancellation arises from the fact that the displacement produced by the long ki 
mode is universal, and hence when the individual hard momenta add up to a soft momentum k, 
the displacement terms also add up as effectively moving a long wavelength mode. Note that the 
cancellation happens independently at any time-slice. Next we show this to hold more generally 
whenever one or several ingoing soft modes are added to a diagram that combines a set of hard 
modes into a soft mode. 

Let us focus again on attaching one soft line via the v ■ V vertex at a specific time r]i, with 

V = (C.5) 

an order m soft velocity, and freeze the rest of the diagram. Cutting this frozen diagram at time 
rji and stripping off the soft ki line we have a set of fields {^^'\qi, f]i)} whose contribution to the 
outgoing field is 

Aaai-{v]Vi,{qi})Yl'^':^^^Hqi,Vi), (C. 6 ) 

i 

where Aa... is a function that contains all momentum dependences arising from the vertices after 
r]i and all propagators. Now add the soft line and sum over all possible attachments of • V vertex 
at rji 

Aaai-iv, m,{qi,--- ,qi + ki, - r?i), (C.7) 

i ^ i 

and use the fact that ^^...( 77 , 771 , {gi, • • • , + fei, • • • } is regular when ki <C Qi to get 

Aaa^-iq, 9 l,{qi,--- ,qi,---})Y If ^l) + (C-8) 

i ^ i 

The sum gives ki ■ {k — ki)/kf = 0 {ki ). Hence, the drifts caused by a long wavelength mode add 
up at any moment to the overall drift of the total (stripped) diagram whose momentum is k — ki. 
This proves that if we start from a composite operator such as which combines hard modes 

into a soft outgoing line, and perturbatively add soft modes they will not introduce any negative 
power of ki- 
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D Equations in Fourier space 


In this appendix we provide explicit formulae for the equations of motion of momentum and velocity 
in Fourier space. 

The equations of motion using momentum are 


dr5 + fj, = 0 , 

+ + + didj 

driy^ + Hv^ = + + 

Several definitions are in order. comes from the term ddicf) 


-u 


mk 


(D.l) 

(D.2) 

(D.3) 




tlm f a{q,k - ^6{^6{k - ^ 

Jq 

( a^^\q,k - ^5{^5{k - ^ , 
Jq 


where a was given in (2.9) and the arrow refers to symmetrization 


12 


= 7 a{k,^ ±a{q,k) 


1 


{k + ^ ( q ^ k 


q 


2 j^2 


(D.4) 


(D.7) 


The other interactions come from di and are as follows. 


e^\k) 

r^Hk) 


E 


n=0-^9l’g2 


^^-^T^2^k‘{q2)p{qi2) - qi){-'^T, 

qku 


„=0 -^91.92 


fc • 92 X 1/(92) k ■ ^L2li(^L2) -k-qu^ i'iqu) k ■ 921^(92) 


9i9?2 


<5-(A:-9-'i)(-1)", 


n=0-^9l>92 


A: • 92 X P{q 2 )k ■ 912 x 1/(912) 


2 2 

qiqi2 


5-{k-q,)i-iy 


E 


^=o■^9l,g2 


fc X g2fc s^{k - 9 i)(-l)” , 

92912 


„=0 •^91,92 


A: X [92 X 1/(92)] A; ■ ^l2/u(9i2) - A; x [912 x 1/(912)] A; • 92 ^( 92 ) 


9i9?2 


r(A:-9-'i)(-l)", 


n=0-^91,92 


A; X [92 X 1/(92)] k X [912 X 1/(912)] 


qhu 


6-{k-qi)i-ir, 


^^The reason why we drop the antisymmetric term is the following. For any anti-symmetric function A{k, q) = 
—A{q, k) one has 

f A{q,k-^ = - f A{q,k-^, (D.5) 

«/ q «/ q 

as it can be easily proven. This equation has two solutions: zero and infinity. In applications to perturbation theory 
the integral is often divergent. We regulate the divergences in a way that respect the anti-symmetry so that the 
regularized terms vanish. E.g. a divergence in a integral like 


f a{q,k - ^5%5{k - fj, 
J a 


(D.6) 


is regulated by smoothing each 5{k) on some cutoff scale fccutoff- 
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where qi 2 = — q 2 , 6^{k) is the Fourier transform of 6{x)^ and we used the vector identity 

a X b X c = ba ■ c — ca ■ b (D.8) 

Moving on to the velocity equations, we start with the Fourier-space decomposition 




(D.9) 


The equations of motion in terms of velocity gradient and curl in Fourier space, supplemented by 
the EFT corrections, are then given by 


dr6 + 9 = / -a{q,k - q)9{q) + aw{q,k - ^ ■ w{q) 5{k - q), 

Jq L -I 


(D.IO) 


drO + ne + 


/3{q, k - ^9{^9{k - ^) + (3w{q, k-^ ■ w{k - ^9{^+ 


w{^r^{q,k- ^w{k- ^ 


+ di 


1 + 5 


drW + l-Lw = f ^^j{q,k - q)w^{k - q)9{q) + d^{q,k - ^w{q) ■ w{k - ^ 

iq 

+7.7(9. k - q)w\k - q)-aP{g) + 
where a and /3 were given in (2.9) and^^ 
aw{qi,q2) = 


(D.ll) 


(D.12) 


q2Xqi 3 qixq2 

— 2 —’ Pw{qi,q2) = - 


qi 


2 ’ 

qi 




(+L X q2)i {qi X q2). 


qhi 


iij{qi,q2) = 


qi,i {qi + q2)j {qi + ^ 2 ) 


qi 


_ d^'\qi + q2)vqi 




Qi 


^ij, ifji{qi,q2)^ 


qI 


(D.14) 


E Glossary 

In this paper, we have extensively used the same nomenclature as in quantum and statistical field 
theory to stress the similarity with the problem at hand and to make the discussion more intuitive. 
However, when needed we have adapted and refined some of the definitions. For the convenience of 
the reader we have collected below those technical terms used in this paper whose meaning differs 
in some important or subtle way from the common usage. 

used the magic vector identity 

d X {v^dk^ = —d X X • iT^ j . (D.13) 


39 
















1-Particle Reducible (IPR) A loop diagram for a correlation function that can be divided 

into two disconnected parts by cutting an internal line that 
is not immediately connected to a contraction or an external 
line, e.g. (3.4) 

1-Particle Irreducible Any loop diagram that is not IPR, e.g. (3.1) 

Time-evolution diagram It shows the evolution of a set of initial fields into one higher 

order final field at a later time (3.6) 

Hard line In a diagram, a line representing the time evolution of per¬ 

turbation 5{k) with a large momentum k ^ oo, e.g. thick 
lines in (3.6) 

Soft line In a diagram, a line representing the time evolution of per¬ 

turbation 6{k) with a small momentum, e.g. thin lines in 
(3.6) 


Amputated diagram 


Non-stochastic 


Stochastic 


Double softness 


A time evolution diagram in which the soft lines are imme¬ 
diately connected to the hard lines, e.g. (3.11) 

The counterterm of a time-evolution diagram in which all 
of hard lines are paired and contracted with each other, e.g. 
(3.8) 

The counterterm of a time-evolution diagram in which one 
or more hard ingoing lines remain unpaired, e.g. (3.9) 

A field is double soft iff it scales as as /c —)• 0, e.g. 6 and 
Vtt 


Systematic Renormalization is systematic if new counterterm are needed 

only for IPI diagrams, while divergences in IPR diagrams 
are already canceled by lower order counterterms 
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